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Abstract 

Molecular dynamics (MD) simulations of crystalline poly(ethylene oxide) (PEO) 
have been carried out in order to study its vibrational properties. The vibrational 
density of states has been calculated using a normal mode analysis (NMA) and also 
through the velocity autocorrelation function of the atoms. Results agree well with 
experimental spectroscopic data. System size effects in the crystalline state, studied 
through a comparison between results for 16 unit cells and that for one unit cell has 
shown important differences in the features below 100 cm~^. Effects of interchain 
interactions are examined by a comparison of the spectra in the condensed state to 
that obtained for an isolated oligomer of ethylene oxide. Calculations of the local 
character of the modes indicate the presence of collective excitations for frequencies 
lower than 100 cm~^, in which around 8 to 12 successive atoms of the polymer 
backbone participate. The backbone twisting of helical chains about their long axes 
is dominant in these low frequency modes. 
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1 Introduction 



Solid polymer electrolytes (SPE), composed of inorganic salts solvated in solid, high 
molecular weight polymer matrices, have been the focus of intense theoretical and 
experimental research because of their applications as solid state batteries [1, 2]. In 
spite of their wide technological applications, the precise mechanism of conduction in 
these materials is still unclear and remains a pursuit of interest. In non-polymeric 
crystalline and glassy electrolytes, the ionic species hop from one site to another 
within a rigid host frame. However, the conduction mechanism in solid polymer 
electrolytes is believed to be different [3, 4, 5]. NMR studies of line shape and 
relaxation rates in poly(ethylene oxide)-lithium salt complexes have demonstrated 
a relationship between the motion of the Li+ ions and the segmental motion of 
the PEO chains [6]. This notion gains support from other experimental studies 
such as Field-Gradient NMR spin echo technique [7], and Quasi-elastic Neutron 
Scattering (QENS) [8]. SPE, in general, contain both amorphous and crystaUine 
regions, with conduction being facile in the amorphous regions. This is probably 
due to differences in the chain dynamics [9]. The work of Armand and coworkers 
has shown that ion mobility in solid electrolytes is of a continuous, diffusive type, in 
the amorphous regions [10]. NMR, differential scanning calorimetry and electrical 
conductivity studies have demonstrated that both cations and anions are mobile in 
the amorphous phase [11]. Recently, Bruce and coworkers have shown that a SPE 
with a ratio of ether oxygen to lithium of 6:1, exhibits higher conductivity in its 
crystalline phase as compared to its amorphous phase and that the ion transport is 
dominated by the cations [12]. In such complexes, cross-linking between a pair of 
polymer chains results in a channel like structure. The ion moves in this channel, 
aided by the segmental motion of the polymer chain pair. 

A large number of simulations have been carried out on PEO and PEO-salt 
complexes, both in their crystalline and amorphous phases over the last decade. 
Molecular dynamics studies by de Leeuw et al, using an united atom model (UAM) 
for PEO chains in their molten state, have examined the motion of methylene groups 
which they term as segmental motion [13]. Neyertz et al have performed extensive 
molecular dynamics simulations of bulk crystaUine PEO [14, 15], PEO melts [16], 
crystalline Nal-PEO [17] and amorphous Nal-PEO [18] systems. The simulations 
of Muller-Plathe and coworkers on PEO-LiI complexes [19] have demonstrated that 
the ether oxygens belonging to consecutive monomers of a PEO chain coordinate 
to the same Li"*" ion and that this segmental coordination is argued to be the driv- 
ing force for salt dissolution. Laasonen and Klein have performed MD simulations 
of both crystalline and amorphous PEO-Nal complexes [20]. Their study showed 
that ion pairing is most probable in the crystalline rather than in the amorphous 
phase. Halley and coworkers have studied the structure of amorphous PEO using the 
Parrinello- Rahman method recently [21]. MD simulations have also been employed 
by Smith and coworkers to elucidate the structure and dynamics of poly(propylene 
oxide) melts [22], PEO-salt complexes [23] and aqueous PEO solutions [24], using a 
potential model derived from ab initio calculations [25]. 
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It is our endeavor here to characterize the vibrational modes of the polymer back- 
bone in the pristine polymer, with an emphasis on the low frequency modes. This 
could enhance our understanding on the exact relationship between the dynamics 
of the polymer and that of the cation in the PEO-salt complexes. We carry out this 
study with the thought that understanding the evolution of the vibrational modes 
from the crystalline phase of pure PEO to its nature in the glassy state without and 
with the salt is crucial in mapping the precise mechanism of ion transport in these 
systems [26, 27]. Normal mode analyses to characterize vibrational spectra have 
been employed successfully to describe the dynamics of glassy systems [28] , and are 
hkely to be employed in the study of jammed granular materials [29]. These calcula- 
tions have also helped our understanding of vibrational modes of polyethylene [30], 
of proteins in solution [31], and in determining the flexibility of proteins, and in the 
thermodynamics of hydration water in protein solutions [32]. Unlike polyethylene, 
PEO adopts a distorted helical conformation in its crystalline state, with the distor- 
tion arising from strong intermolecular interactions. It is thus important that such 
normal mode calculations be performed for the crystalline state, rather than for an 
isolated oligomer. In this paper, we limit ourselves to the study of the vibrational 
spectrum of PEO in its crystalline state. Analyses of these modes under other state 
and phase conditions will be examined in future. In anticipation of our results, 
we have identified the existence of segments in the polymer backbone consisting of 
around 8 to 12 atoms, to exhibit significant atomic displacements, for vibrational 
modes with frequencies up to around 100 cm~^, and that the modes with frequen- 
cies below 60 cm^^ involve twisting of the skeletal backbone. The paper is divided 
as follows. Details of the simulation and the methods of analyses are presented in 
the next section. Later, we discuss the results obtained from our study. Details on 
obtaining the elements of the dynamical matrix are provided in the Appendix. 



2 Details of Simulation 

The unit cell of crystalline PEO is monoclinic with the PEO chains in a (7/2) dis- 
torted helical conformation with TTG sequence [33]. An unitcell consists of four 
PEO chains with 21 backbone atoms in each chain. The MD runs reported here 
were initiated from configurations generated from this crystal structure. The two 
ends of a polymer chain were assumed to be bonded across the simulation cell bound- 
ary [14, 20] to eliminate end effects. MD simulations were performed primarily for 
a system that we describe as (422), which contained 4 unitcells along the a-axis, 
2 unitcells along the b-axis and 2 unitcells along the c-axis, containing a total of 
3136 atoms. To study system size effects and effects of inter-chain interactions, we 
have performed additional simulations of only one unitcell, and also of one isolated 
finite length polymer chain. Thus the simulations for the (422) cells and that for the 
one unitcell contained chains that had no ends, while that for the isolated molecule, 
contained a chain with two ends. An all atom model (AAM), with explicit consid- 
eration of hydrogen atoms was used, to properly account for the steric interactions 
expected to be dominant in the crystalline state. The simulations were carried out 
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in the canonical ensemble at 5K for the NMA and in the constant pressure ensemble 
at 300K and 1 atmosphere to test the stability of the simulated crystal. Tempera- 
ture control was achieved by the use of Nose- Hoover chain thermostats [34], using 
the PINY-MD program [35] . Long range interactions were treated using the Ewald 
method with an a value of 0.3 A~^, and 2399 reciprocal space points were included 
in the Ewald sum [36]. The methylene groups were treated as rigid entities, en- 
abling us to employ a timestep of 1 fs. To obtain the vibrational density of states, 
we performed these calculations at a temperature of 5K, where the atoms could be 
expected to be near their equilibrium positions. The equilibration period for the 
single molecule, the unitcell, and the (422) system were 65 ps, 100 ps, and 750 ps, 
respectively. These were followed by an analysis run of duration 30 ps during which 
the coordinates and velocities of each particle were stored at regular intervals. Ve- 
locities were dumped every time step to obtain the power spectrum of their time 
correlation functions and the coordinates were dumped every 10 fs. 

The simulations were performed with a force field obtained from the work of 
Neyertz [14] with the torsional parameters of the CHARMm model [37]. The poten- 
tial parameters are given in Table 1. The non-bonded interactions were truncated 
at I2A for the (422) crystal, at 3.25A for one unitcell, and at I2A for the single 
molecule runs. The lower cutoff value for the unitcell run is necessitated by the fact 
that the interaction cutoff should be less than half of the minimum distance between 
any two opposite faces of the simulation cell. 

The potential energy of the system can be expanded in terms of atomic displace- 
ments from an equilibrium configuration as, 

^ rr/ ^ (dU\ I ( d^U \ 

C/(gi,g2,-,gn) = f^(goi,go2,-,gon) + ^ ^» + o a o + - (1) 



%/q 2 \dqidqj 







where the subscript represents the equilibrium configuration and rji are the devi- 
ations of the coordinates from equilibrium, represented as, 

Qi = qoi + rji (2) 
The elements of the Hessian matrix are given by 

'•^ y/rriimj \dj3jdai j 

where i,j represent particle indices and a, {3 represent the spatial coordinates x,y,z. 
rrii is the mass of particle i. A simple scheme to obtain some of these Hessian 
elements efficiently is provided in the Appendix. All Hessian elements obtained from 
such analytical expressions were checked against numerical second derivatives [38] 
within our normal mode analysis code, and were found to match. The eigenvalues 
and eigenvectors of the Hessian matrix were examined to understand the vibrational 
dynamics of PEO. The frequency, Ug ,of a particular mode of vibration, s, is related 
to its eigen value, , by 

A. = (27ri.,)' (4) 
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With these set of interactions, the initial pressure of the (422) system was found 
to be around 4000 atmospheres. Hence, we performed a MD run in the NPT en- 
semble for 200 ps under ambient conditions. The change in volume was found to be 
1.3% from that of the experimental crystal. Time correlation functions of atomic ve- 
locities were calculated, and were Fourier transformed to obtain the power spectra. 
These were compared with the spectrum obtained from the NMA. The spectra were 
convoluted with a Gaussian function of width 4cm~^ so as to provide a width to 
the spectral features. The helical axis of a PEO chain possesses a 7-fold rotational 
symmetry. In addition, 7 C2 axes lie perpendicular to it, making the molecular sym- 
metry of PEO to be D7 [33]. We have characterized the symmetry of the normal 
modes by studying the transformation of atomic displacements, on application of 
the symmetry operations for this group. Specifically, the operation by the C2 axes 
enables a distinction between the three irreducible representations, Ai, A2, and E. 

We have also characterized the spatial extent of the modes of vibration using a 
quantity called the local character, defined as [31, 39] 



where Uij is the i component of the f eigenvector. The local character value can 
range from to 1 and it determines the extent of localization of a particular mode. 

To obtain quantitative information on the length of the segment involved in the 

low frequency modes, we have defined a quantity called the Continuous Segment 
Size (CSS). We calculate the displacement of the k^'^ atom due to a i*'* mode using 
the expression 



where ruk is mass of k^^ atom, {Tf is the vector formed by the components of the i*'* 
eigen vector contributed by the k*'* atom, T is temperature, and cui is the frequency 
of the i^^ normal mode. We then check if the displacement is greater than a specified 
cutoff. If n successive backbone atoms of a chain each have displacements greater 
than the displacement cutoff, they are defined to constitute a segment with CSS = 
n. Similarly, CSS was calculated for all possible modes with different vibrational 
frequencies from which the average segment size for a given frequency was calculated. 

3 Results and Discussion 

For crystalline systems, a close match between the experimentally determined cell 
parameters and that obtained from simulations is a crucial first step in the veracity 
of the parameters used. We show in Figure 1, the time evolution of the cell param- 
eters at 300K and 1 atmospheres, generated by a MD run in the constant pressure 
ensemble using the Parrinello- Rahman method [40]. The b and the (3 parameters 
of the unit cell exhibit a relaxation from the zero time experimental value, due 



3N 

Localcharacter (j) = Xl^t? 

i=l 



(5) 



< 5rl kBT 




(6) 
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possibly to relatively minor deficiencies in the interaction model used to represent 
this system. The constancy of the cell parameters for over 100 ps, shows that the 
simulated crystal is in a stable state, and that the potential parameters are indeed 
able to reproduce the high demands of the crystal symmetry. Table 2 compares the 
cell parameters obtained from our simulations to the experimental data. 

An interesting feature associated with polymer dynamics is the variation of the 
various vibrational modes of a chain encountered during its transformation from the 
isolated state to the crystalline state. The vibrational spectra of a single molecule, 
one unitcell and that of the (422) crystal are compared in Figure 2a. When the 
chains assemble to form a crystal, each chain might prefer a new conformational 
state relative to its structure in the isolated state. A comparison of the vibrational 
spectra between that of one isolated molecule of finite length and of the (422) sys- 
tem provides information on the effect of intermolecular interactions. As expected, 
the vibrational states of the isolated molecule showed marked differences from the 
crystalline state, particularly in the low frequency regions, where large amplitude, 
collective motions are predominant (see later). However, there are no significant 
changes in the high frequency regions of the spectrum. The spectrum for the one 
unitcell compares well with that of the (422) crystal. However, the features in the 
VDOS of the (422) crystal is much better resolved, particularly at low frequencies. 
For instance, the features at around 40 cm^^ and 75 cm^^ are clearly evident in 
the larger system than in the spectrum for the unitcell, pointing to effects of long 
range interactions. The spectrum for the (422) system is shown in an expanded 
scale in Figure 2b. The split in the feature below 100 cm^^ is evident and compares 
well with experimental infrared (IR) spectra [41] that shows features at 37 cm~^, 
52 cm~^, 81 cm"^ and 107 cm~^. We do observe a prominent shoulder at 108 cm~^ 
which has been attributed to modes involving C-0 internal rotation earlier [41]. A 
comparison of the simulated vibrational density of states with the experimental IR 
and Raman spectra, exhibited in Figure 2c, shows that the potential model captures 
well nearly all the vibrational modes [41, 42, 43]. Note that the simulated spectrum 
is the raw density of states and does not contain any other terms that are needed to 
calculate the experimental spectra. Thus only the peak positions need to be com- 
pared, and not their intensities. Almost all the features found in the IR and Raman 
spectra seem to be present in the simulated VDOS. Another noteworthy feature of 
the spectrum is the absence of modes with imaginary frequencies which would have 
corresponded to either the presence of atoms away from equilibrium locations or to 
a mismatch between the empirical potential function and the crystal structure. 

We have also calculated the vibrational spectrum through the Fourier transfor- 
mation of the velocity auto correlation function of the atoms. This is compared 
with the VDOS obtained by the NMA in Figure 3. The agreement between the 
spectra is excellent except for features above 1200 cm^^. Also, the peak at around 
1470 cm~^ is missing in the spectrum obtained through the velocity autocorrelation 
function (VACF). Visualization of the atomic displacements associated with this 
mode revealed HCH bending in methylene groups. Since all the CHg groups were 
constrained to be rigid during the MD runs for the VACF analysis, the absence of 
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a peak in the VDOS obtained from VACF, in this region can be rationahzcd. The 
comparison of the VDOS obtained from the two methods is good, despite the fact 
that the NMA method is only a harmonic approximation to the potentiaL Close 
examination of the two shows a marginal (approximately 4-5 cm"-*^) shift to higher 
frequencies in the spectrum obtained by the NMA method relative to that from the 



We visualized the eigenvectors corresponding to different vibrational modes to 
assign the nature of atomic displacements that are responsible for the spectral fea- 
tures. In general, our assignments are consistent with earlier calculations [44, 45]. 
In Figure 4a and Figure 4b, we display a few of the modes. The zero frequency 
mode characterizes rigid body translation. The features at around 40 cm~^, have 
earlier been attributed to chain deformations [41]. Based on visualization of atomic 
displacements shown in Figure 4a, we assign these modes specifically to the twist- 
ing of the polymer backbone. The mode at 88 cm~^ arises from torsional motion 
around the C-0 bond. The Ai mode at 216 cm~^ can be assigned to torsions around 
the C-C bond, while the 510 cm~^ feature involves bending of CCO triplets. COC 
bending is observed at 952 cm~^ while the wagging motion of the methylene groups 
is found to be present at 1244 cm~^, in good agreement with IR and Raman mea- 
surements [44]. The symmetry of the modes are also shown in Figures 4 agree well 
with experimental assignments [44, 41, 45]. 

For characterizing the normal modes further, we have calculated the local char- 
acter for each mode [31, 39] which is shown in Figure 5. In the range, to 160 cm~^, 
the local character value is very small, implying the participation of a large number 
of atoms. However, the local character does not provide any information on the 
proximity of the atoms that exhibit significant displacement in a mode. Thus this 
quantity has to be augmented by a further analysis of the concomitancy or proximity 
of atoms excited in a mode. 

We have developed such an index that takes into account the connectivity of the 
atoms involved in a mode. We determine the number of successive atoms, n, that 
participate in a mode of a given frequency, by calculating the distribution of segment 
sizes, f(n), for that vibrational mode. As a representative example, the distribution 
of the CSS for one such vibrational mode of frequency 44 cm^^ is shown in Figure 6 
for four selected displacement cutoffs. It is our contention that a large number of 
atoms in proximity to each other participate in these modes. This is evident from 
the non-zero value of f(n) for n values in the range of 5 to 10, for some of the cutoff 
values. 

The average segment size was calculated by evaluating the following summation. 



In our study, we define a segment as a chain of connected atoms whose displacements 
are larger than a specified cutoff, with n > 5. Values with n < 5, which correspond 
to intramolecular excitations that arise out of bonded interactions like the stretch, 
bend and torsion, have been omitted in < CSS > calculations, as their origin is 



VACF. 



< CSS > 




(7) 
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trivially known. It should also be noted that the exact values of the average segment 
size will depend on the chosen displacement cutoff. Hence, we have performed these 
calculations for a variety of cutoffs and these are exhibited in Figure 7 as a function 
of the frequency of the modes. Notice that the segment size for the zero frequency 
modes, i.e., translations, is 42, which is the number of backbone atoms in a given 
chain. It is also evident from the figure that the average segment size decays with 
increase in frequency. Figure 7 also shows the variation of <CSS> as a function 
of frequency, for four different displacement cutoffs. The displacement cutoff which 
spans the frequency range relevant to collective motion, where the the local character 
value is almost zero is the one of significance. It can be seen that in the frequency 
range of 10 cm"^ and 100 cm~^, which is the range of collective motion, around 8 
to 12 successive atoms of the backbone are involved in the excitations. 

4 Conclusions 

We have studied the vibrational dynamics of crystalline poly(ethylene oxide) using 
molecular dynamics and normal mode analysis. The vibrational density of states 
obtained from NMA matches well with that obtained from the Fourier transforma- 
tion of velocity autocorrelation function and also with experimental IR and Raman 
data [41, 42, 43]. The VDOS of an isolated PEO chain of finite length showed 
marked differences from that of the crystal in the low frequency region where collec- 
tive modes are predominant. We have also explored system size effects by comparing 
the VDOS obtained from a simulation containing 16 unitcells and that of one unit- 
cell. The spectrum obtained from the former is much better resolved, particularly 
at low frequencies, where three clear features, at 44 cm~^, 70 cm~^, and 88 cm~^ 
are observed. The results of our calculations agree well with assignments of mode 
symmetry by earlier workers, that used standard methods for a single chain of PEO, 
or for an unit cell [44, 41, 45]. Such calculations have the added advantage, over the 
MD route presented here, of being able to obtain dispersion of the vibrational modes. 
However, the method described here can be used to characterize these vibrations in 
the amorphous and liquid phases of PEO. 

The normal modes obtained from the present analysis were characterized by 
the local character indicator and by a new quantity that determines the number of 
concomitant atoms excited by a mode, called the CSS. In the range to 160 cm~^, 
the value of the local character indicator was very small, indicating the participation 
of a large number of atoms in the vibrational modes in this range. A distribution 
of segment sizes was calculated for each mode from which the average continuous 
segment size was calculated as a function of the vibrational frequency. The frequency 
dependence of <CSS> clearly shows collective modes to be present for frequencies 
less than 100 cm~^, in which around 8 to 12 successive atoms of the backbone 
participate. This quantitative analysis is also corroborated by visualization of the 
atomic displacements for the low frequency modes. 

The phase behavior of PEO and the evolution of these vibrational modes as a 
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function of temperature and their properties in the amorphous phase will form the 
objectives of our study in future. 
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6 Appendix 1 



The form of the potential used in our simulations is generic to macromolecular 
interactions, and can be found in Ref . [14] . Here, we provide only the torsional and 
Coulomb terms. 

jj _ QiQj /o\ 

'^coulomb — y'-' ) 

Tij 

and 

6 

Utorsion = f {cOS(p) = ^ OjCOsV (9) 
i=0 

We provide here a procedure to obtain the Hessian elements, which is easy to pro- 
gram. To our knowledge, such a scheme has not been outlined so far [46], and hence 
we provide it here for pedantic reasons. We limit these details to contributions 
from the torsional interactions, the reciprocal space part and the real space part of 
the Ewald sum. Contributions from other terms in the potential energy are much 
simpler to derive and are not given here. 

6.1 Hessian from the torsional interaction 

The torsional angle (p between the planes formed by the bond vectors fi2, f23, and 

^34 is, 

coscf) = i • S = A^B^ (10) 
where the Einstein summation convention has been used and 

A=-^ f^, and B = ^ (11 

I ri2 X r23 I I r23 x r34 I 

The second derivative of Utorsion with respect to a spatial coordinate «„, where a 
can be either x, y, or z and n=l, 2, N, can be written as. 



d"^ Utorsion _ ( df \ f d'^coscp \ ^ f (9cog0 \ / dcos(p \ f f 

OPrndan \dcOS(f) ) \d/3mdan J \ dan J V 9Prn J \dcOS(j)'^ 



(12) 



The derivates of cosd can be calculated as follows. 



d'^cos(J) . d'^B^ dB^dA^ ^ d'^A^ dA^dB^ . , 



dPrndari d(3mdan dari dPrn dfimdan dan dPrr 



We can write A^ as 



with 

tA 



K = «7«£ (ft^). (n!3)j (15) 

and 

2 



Da = 



(16) 
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where 7, v, and ^ stand for any of the three indices x, y, z and e^^^ is the anti- 
symmetric Levi-Civita tensor of rank 3. A similar expression can be written for 



By. 



The first derivative of are, 



1 aTV^ 

'a 



dNf 



where, 



dar^ ^ 



dan Da dan ^ ' dan 

iv<x Wvi)^ [bnZ - Sn2) - {r23)^ (<^n2 " ^nl)] 



The second derivative of is calculated as follows. 



d dA. 



d 



■7 



1 dN:^ 



Da dan 



d 



dPrn dar. 



dPn 



dPrr 



with 



d 



1 



Da dctn 



1 d dN^ 1 dN^OD, 



d 



dPrr 



m dar. 



DAdPrn dan D\ dan dPn 



'a 



dPn 



7 



E 



Nr 



'dPrndar 



+ 



'dN,\ ( dNi 



da. 



m 

dan 



' 1 dNy 



dPm^ 

dDA 



TV 



where. 



d dN"^ 

TTZ n ~ ^7/8a! [{Sm2 — ^ml) (<^n3 " <Jn2) — (^m3 — ^^2) (^n2 — (^nl) 

dpm dan 



(17) 



(18) 



(19) 



(20) 



(21) 



(22) 



Using equations (13), (17), and (19) we can calculate which can be substi- 

tuted in equation (12) to get the torsional contribution to the Hessian. 



6.2 Hessian from the Coulomb interaction 

The reciprocal space energy in the Ewald sum method is. 



where. 



ak — ^2 ?jCos (Jc ■ fj^ , and &fe = E ' 



(23) 



(24) 
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Here, V denotes the volume of the simulation cell, ( determines the width of 
the Gaussian charge distribution centered on the point charges in the Ewald sum 
method, k denote the reciprocal lattice vectors, and eo, the permittivity of free space. 

It can be shown that 



d^Ur 



dPrndar, 



— K 

— -p- [^qnka -kpSmn {cikCOS (k ■ f„) + bksin 



+ QrukpCOS (k ■ ifn - Tm))] } 

The expression for the real space energy in the Ewald summation is, 

Ureal 

where. 



47re, 



N N 



i=i j>j 



er/c (Cnj) 



The second derivatives of U^eai are 

d ( dUreal \ 



N N 



d f dp 



with 



dPrn V dan J 47reo ^ 
d f dp \ [ dp \ d I drij\ 



+ 



dvij \ d ( dp 



dPrn \danj \drij ) dPrn \doin I \dan ) dPrn \dr, 



d ( dVij \ {Sjn - Sin) {Sjm " 4 



It can be shown that 



Sa3 — 



dp_ 
dUj 



= - (si + S2) 



and 



with. 



2 (si + S2) 



Si 



2^6^^ 



and S2 



erfc (Cnj) 



These can be used in Eq. 28 to obtain the Hessian elements. 



(25) 



(26) 
(27) 

(28) 

(29) 
(30) 

(31) 

(32) 

(33) 
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Table 1: Parameters of the interaction potential [14, 37]. 



h [K A- 



Stretch 



ro 



C-C 

C-0 
C-H 
H-H 



1.53 

1.43 
1.09 
1.78 



237017.0 

171094.8 
Constrained 
Constrained 



Bend 




ke [K rad-2 




C-O-C 




112" 






110255.50 






0-C-C 
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0-C-H 




109.5" 






30193.20 
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45199.50 
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ao [K] 


ai [K] 


02 [K] 


as [K] 


a4 [K] 
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as [K] 


C-C-O-C 


-50.322 
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0.0 


-201.288 
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C-O-C-H 
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Atomic charges 
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Table 2: Lattice parameters obtained from simulation compared to experiment [33]. 



Lattice parameter 
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Experiment 


a[l] 
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c[A] 
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125.40 


in 
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Figure lb: 

Figure 1: Plot of instantaneous cell parameters of PEO crystal obtained from sim- 
ulations : (a) cell lengths, (b) cell angles. The zero time configuration corresponds 
to the experimental crystal structure. 
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Figure 2 a: 



20 




600 



V (cm ^) 

Figure 2b: 



21 




Figure 2c: 

Figure 2: (a) The vibrational density of states of a single molecule (bottom), for 
one unitcell (middle) and for the (422) crystal (top) of PEO, each obtained from 
normal mode analysis. The two sections of each spectrum are normalized indepen- 
dently to enable better comparison among the three system sizes. The spectra in 
Figures 2a, 2b, and 2c are convoluted with a Gaussian function of width 4cm~'^. (b) 
The vibrational density of states of the (422) crystal, obtained by NMA shown in 
expanded scale, (c) Vibrational density of states obtained from the NMA method 
(bottom panels) are compared with experimental Raman (Dashed lines) and Infra- 
red (Continuous lines) absorption intensities. Experimental data presented in the 
top left and top right panels were obtained from Ref. [41] and Ref. [44] respectively. 
Note that the simulated spectrum is the raw density of states and docs not contain 
any other terms that are needed to calculate the experimental absorption spectra. 
Thus only the peak positions are comparable, and not their intensities. 
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Figure 3: The vibrational density of states of (422) PEO crystal obtained from nor- 
mal mode analysis (top panels) is compared with that obtained as the power spec- 
trum of the velocity autocorrelation function of all atoms (bottom panels) . The two 
sections of each spectrum are normalized independently to enable better comparison 
among the two methods. The spectra are convoluted with a Gaussian function of 
width 4cm~^. 
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Figure 4a: 
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Figure 4b: 



Figure 4: (a),(b) : Atomic displacements of representative normal modes of the 
PEO crystal. For each mode, one of the chains which exhibit significant atomic 
displacement is shown, for clarity. Red (or Black) spheres denote carbon atoms and 
white spheres denote oxygens. Hydrogen atoms are not shown for clarity in all the 
modes, except the one with frequency 1244 cm~^. Arrows denote the directions of 
atomic displacements and their lengths are scaled up for the purpose of visualization. 
The mode at cm~^ is the rigid body translation of the chain, and that at 1244 cm~^ 
arises from wagging of CH2 groups. Frequencies in cm~^ are as provided in the 
figures. The representation of these modes in the D7 group, are as follows, with 
frequencies in cm~^ given in paranthesis: A2 (38), A2 (44), E (88), Ai (216), E 
(510), E (952), and E (1244). 
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Figure 5: Local character indicator for the normal modes of PEO crystal. The inset 
shows the frequency range where collective motion is most probable with very small 
local character indicator values. 
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Figure 6: The segment size distribution, f(n), corresponding to the mode at 44 cm~^ 
of the PEO crystal. The distribution is shown for four different displacement cut- 
offs. Such distributions are used to calculate the average continuous segment size 
(<CSS>) which is defined in equation (7). 
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Figure 7: Variation of the average continuous segment size (<CSS>) with respect 
to the vibrational frequencies for different displacement cutoffs. 
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